Most singular vortex structures in fully developed turbulence 
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Using high Reynolds number experimental data, we search for most dissipative, most 
intense structures. These structures possess a scaling predicted by log-Poisson model for the 
dissipation field s r . The probability distribution function for the exponents a, e r ~ e aa , has 
been constructed, and compared with Poisson distribution. These new experimental data 
suggest that the most intense structures have co-dimension less than 2. The log-Poisson 
statistics is compared with log-binomial which follows from the random /3-modcl. 

PACS number(s): 47.27. Ak, 47.27. Jv 



Self-similar properties of turbulence, suggested by Kol- 
mogorov [1], have been intensively studied for a long 
time. The theory predicted simple scaling for the lon- 
gitudinal velocity increments, (\u(x + r) — u(x)\ p ) ~ r^ p , 
where £ p = p/3. It became clear, however, that there are 
corrections to this scaling, £ p = p/3 + r(p/3), due to in- 
tcrmittency. A theory, incorporating the intermittency, 
the refined similarity hypothesis [2] , links the statistic of 
these corrections with the statistic of the dissipation field 
E r , the energy dissipation averaged over a ball of size r. 
Namely, (ef) ~ r T ( p \ Many models have been proposed 
to explain intermittency. It was originally suggested that 
the statistics of e r is log-normal [2]. More recently, She 
and Leveque [3] (hereafter SL, see also [4], [5], and re- 
cent study [6]) have proposed log-Poisson statistics for 
the dissipation field, the theory resulting in a remarkable 
agreement with the experimentally found £ p in [7], [8]. 
These experimental exponents are obtained in Extended 
Self-Similarity approach, which is useful because of ex- 
tended scaling range. In the following, we will also use 
these experimental exponents to compare with SL theory 
with slightly modified parameters. 

Consider the Poisson distribution (see, e.g., [9]): 
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and let a < 0, and therefore b = lnmaxe r . Calculating 
the moments, (£?), we note first that, as 



= 1, we get, 

b = £(1 — e a ). Second, by definition, (e£) ~ r T ^ p \ and 
therefore, £ = Cln(^/r), £ being external scale, and C is 
a constant. Finally, denoting 7 = C(l — e a ), we obtain, 



T{p)=C[l-{l- 1 /Cf}-jr 1 . 



(2) 



SL is recovered from (2) if C — 2, 7 = 2/3, so that a = 
In (2/3), C being the co-dimension of most dissipative 
structures, and 7 is defined by the dissipation rate, i.e., 
inverse time-scale, l/t r ~ r ~ 2 / 3 [3]. 

The meaning of C becomes even more clear directly 
from (1): the most intense fluctuations correspond to 
a = (as a < 0), so that the probability P(a — 0) = 
= (r/£) c — (r / 1)( d ~ Hq \ D is dimension of space 



(= 3). Thus the Hausdorff dimension for most dissipative 
structures in SL theory, Hq — 1, i.e., the structures are 
filaments. On the other hand, using expressions for b, £ 
and 7, we now rewrite (1) as follows, 
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Putting a = in (3), we can see that the most intense 
structures are expected to scale ~ r -7 . This scaling pre- 
dicted by the log-Poisson statistic is proved to be possible 
to verify experimentally. 




FIG. 1. Scaling for most intense structures. The power law 
fitting of the experimental data (solid thick line) has been 
extended to reach unity (solid line), where it is supposed to 
match with SL scaling. The distances are given in terms of 
Kolmogorov micro-scale n. 

We used 10 million points of atmospheric data from 
Yale University, with an estimated Taylor microscale 
Reynolds number 9540 (courtesy of Sreenivasan) . The 
data are treated in spirit of Taylor hypothesis, that is, 
the time series is treated as one-dimensional cut of the 
process. We denote u> = d x u, u is longitudinal velocity 
along the x-axis. Therefore, the dissipation, e ~ u> 2 , and 
we will deal with dimensionless dissipation e = uj 2 /{us 2 ). 
In our case, 
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x+r/2 e{x')dx' 



-r/2 



(4) 



and maxima of e r can be measured. Figure 1 shows a 
remarkable scaling for this quantity for 4.5 decades. The 
deviation from SL is small, and we recall that SL suggest 
that there is no anomalous scaling for t r . This small 
deviation in Fig. 1 can be interpreted 
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stead of making use of inverse Legendre transform of the 
spectrum (2), we provide direct measurements of the ex- 
ponents a. 

Figure 2 shows two typical PDF's for two distances. It 
is clearly seen that the SL theory is well confirmed: the 
experimental PDF's are quite close to the correspond- 
ing Poisson distributions. The deviation is observed only 
at large a, i.e., (exponentially) small e r , for which the 
theory does not claim to account for. The latter 
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FIG. 2. Constructing PDF's for experimental and Gaussian processes, and comparing them with corresponding Poisson 
processes. 



as anomalous persistence of the eddies, which is indeed 
observed [10], see also discussion in [11]. The value of 
7 is 0.61 ± 0.01, only slightly smaller than 2/3. Note 
that maxima are not always predictable; e.g., for log- 
normal distribution, the maximum could be anything. 
Another point: the measurements of maxima is mean- 
ingful because the coarse-grain averaging (4) is already 
present. In order to compare with a "regular" ran- 
dom process we generated a Gaussian process u> g with 
correlation function coinciding with experimental, i.e., 
(u) g (x + r)ujg(x)) — (u(x + r)ui(x)). Then, the "dissipa- 
tion" e (f) = uj 2 g , and ei 9) = l/r f*^/2 e (») (x')dx' . Cor- 
responding calculation for the maxima are reported in 
Fig. 1. If any scaling can be extracted from the Gaus- 
sian process, it would be at large asymptotic distances, 
and the scaling is trivial, 7 = 0, meaning no singularity. 

Another testing of the theory is presented by direct 
measurements of the exponents a in (3). We may con- 
struct a PDF for a-distribution, measuring lne r , and tak- 
ing a = ln{2/3}, say. In fact, the calculated value of a is 
only slightly different from In {2/3}, see below, and, as a 
result, the PDF plots with "true" a are actually indistin- 
guishable from these with a = In {2/3}. In order to com- 
pare the experimental PDF with the Poisson distribution 
we find £ from (3), £ = ln{maxe r }/(l — e a ), and hence 
the distribution (1) is unambiguously defined. Thus, in- 



is constructed to account for asymptotically high mo- 
ments, i.e., high values of e r . As to the Gaussian pro- 
cess, it is clear that first, it is dramatically different from 
the experimental, and second, it is quite different from 
corresponding Poisson distribution (which is defined by 

£ 9 = ln{max£r^}/(l — e a )). The latter predicts much 
higher level of strong fluctuations (corresponding to low 
a) than actual Gaussian - which is hardly surprising. 

Returning to the Legendre transform, we note that, for 
comparison, we have to express a! in (1) through Stirling 
formula, to get, 

1 „ Q -g + Q l nf£ /q> _ 1 Q D ~" , (5) 
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where 

D-H = C(l-y+ylny), y = (h+ 7 )/(C\ ln{l - 7 /C}|), 
and h = —7 — a\ ln{l — 7 /C}|/ln{r/^}. In spite of the 
fact that Stirling formula is valid only for a ^> 1, for- 
mula (5) for our parameters (in particular, £ = 14 16) 
is quite close to (1) except for a < 1 (and, according 
to (5), P — > 00 when a — ► 0). Moreover, the exponent 
D — H exactly coincides with that obtained from Legen- 
dre transform, i.e., D — H = D — inf p {ph + D — r(p)}. 
Now, min h = —7, corresponds to maximal excitation, 
according to (3), while the corresponding Hausdorff co- 
dimension D — H(minh) — C [3]. 
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Although the experimental 7 is not that different from 
2/3, the value of the other parameter C is quite sensi- 
tive to that difference. In order to find C we substitute 
7 from our measurements into (2), and use computer 
routines to find a best fit for these data with free pa- 
rameter C and the exponents C,p ESS ^ from experiment 
[7], [8]. As a result, we find C = 1.67 and a = -0.45 (cf. 
In {2/3} = —0.41). With these parameters, the deviation 
of these calculated exponents (^f 1 from the experimen- 
tal exponents, y] {{ffl - C^ SS) ) 2 ) = 0.0063. To com- 
pare: for SL, ^ ((d SL) - C ( p ESS) ) 2 ) = 0.0078. The ( { p e) 

exponents seem to be "better" than (p SL ^ , but consider- 
ing that the experimental exponents have errors of about 
±1% [8], we conclude that these exponents are similar. 
Note that if we substitute in (2) the value of 7 from the 
experiment, and put C = 2, then, for the obtained expo- 
nents, Cp C=2) , we have, \J {{C,¥ =2) - ( { P ESS) ) 2 ) = 0.050, 
much too high. Figure 3 shows r(p/3) from experiment, 
and for different theories. It can be seen that all the 
curves collapse into one, corresponding to the experi- 
ment, except that one with 7 from our measurements, 
and C = 2. This illustrates that the data are indeed sen- 
sitive to the measured 7, that is to its (small) difference 
from 2/3. The codimension C = 1.67 corresponds to 
H n = 1.33. This value of H > 1 seems to be consistent 
with the distinction between persistent vortical filaments 
and the dissipative structures associated with regions of 
strong strain [13]. That means that the most dissipative 
structures consist not only of filaments, but in part of 
sheets, or filaments convoluted into complex structures, 
covering more than 1 dimension. 
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FIG. 3. Intermittency corrections from experiment [8], and 
from other theoretical models. 

Note that, according to the intersection theorem [14], 
dimension H n < 2 cannot be detected in ID measure- 
ments directly, and therefore our conclusion is inevitably 
indirect. It would be important to measure the Hausdorff 
dimension in 3D simulations directly. Another reason for 



that is the surrogacy issue [15]. 

These statements about the dimensions of the most 
intense structures can be also formulated for log-binomial 
distribution, and, it is known that the Poisson process is 
a limit of the binomial distribution for "rare events" . In 
particular, the Poisson distribution can be obtained from 
the random 0- model [12] by a suitable limiting process 
[4], [16]. Let take two values, W — 0\ with probability 
x, and W — 02 with probability 1 — x, and 0\x + #2(1 — 
x) = 1 (in order to have r(l) = 0). Let also x < 1/2, and 
01 < 1 < 02- Then, on the n-th level, the distribution is 
binomial, that is, W n — s n — 0™0% ~ m with probability 
(™ )x m (l - x) n - m . Hence, (e£) = [x0{ + (1 - x)0%] n . 
Taking into account that n — ln(r/£)/lnT, T being the 
ratio of successive scales, we obtain, 

_ In [x0l + (1 - x)0l] 
lnT 

In [4] and [16], T was treated as a free parameter. It was 
shown that, if T = 1 — x/C, 0\ = 1 — 7/C and x — > 0, 
then 02 ~ 1 + x-f/C, and (6) reduces to (2). 
The most intense structures on n-th level, 

/ r \ ln/3 2 /lnr /r\-f/3 

=G) < (7) 

cf. (3). On the other hand, the probability of these 
maxima, 

ln(l-x)/lnr /r\ c & 

p=(i-*r=( 7 ) =y . (8) 

In particular, if T = 1 — x/2 and x — > 0, then Cp = 2 
[17]. 

If we do not treat T as a free parameter, we may con- 
sider T = 1/2 1 /- , say. As 1 < 2 < 2, 7/3 , accord- 
ing to (7), satisfies < 7/3 < D. And, according to 
(8), < C (i (= ln(l -x)/lnT) < D. Thus, both 7/3 
and Cp satisfy requirements for codimensions. A par- 
ticular case x = 1/2 corresponds to the model proposed 
in [18]. Then, according to (8), P = {r/£) D , i.e., the 
Hausdorff dimension H is = 0, while 7^, defined from 
(7), = Dln/3 2 /m2. The case 0\ = returns us to 
the /3-model [19]. In that case, 02 = 1/(1 — x), and 
1p = Cp = D\n{l-x)/ In (1/2). 

The log-binomial distribution generally cannot be re- 
duced to the log-Poisson PDF: in particular, even if x 
is small, and T = 1/2 1 / 13 , then, r(p) ~ x — > 0, and 
thus the intermittency is negligible . Then, at first 
sight, the log-binomial distribution is not suitable for 
our purpose. Indeed, according to the central limit the- 
orem, the binomial distribution is asymptotically nor- 
mal at n > 1, and the log-normal distribution has un- 
surmountable shortcomings [20], [21]. Nevertheless, the 
spectrum (6) does not even look like log-normal (for 
which t(p) = — (u/2)p(p— 1)) and rather behaves like log- 
Poisson for p ^> 1. Indeed, according to (6), for p ^> 1, 
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t( P ) = C /3 + Ci^-p7/j, C-i = 



(1-jc) lnT 



<0, (9) 



/3 = Pi/p2- This spectrum resembles (2); and the con- 
stants in (9) happen to be numerically close to corre- 
sponding numbers in (2). The reason for such a dramatic 
difference with log-normal distribution is as follows. For 
binomial distribution, 



(O=/£ P £(m)s m (l-a0 n_,n 

m=0 



(10) 



S = [3 P <C 1 for large p. Then <5 m decreases dramati- 
cally with increasing m, and therefore the terms of the 
sum (10) of maximal probability, at m ~ xn, where nor- 
mal distribution if formed, do not contribute substan- 
tially. In contrast, only the first few terms of this sum 
(responsible for "rare" and very intense events) really 
contribute. Thus, effectively, the distribution works like 
a Poisson distribution. To see this explicitly, consider 
a probability distribution (™ )x m (l - x) n ~ m 8™ /A, where 
A is a normalization constant, A = (x5 + 1 — x) n , and 
So = {P1/P2Y , Po ^ 1- Then, for large n we express 
the factorials entering the binomial coefficients through 
Stirling formula (except for ml, because m is not neces- 
sarily large), to get, 



Pi(m) = jUx m (l-, } 



n — m cm 



ml 



where £0 = nS n x/(l — x). For p > p , the sum (10) can 
be written as Aft* Em=o P i (m)S m( - p - p °\ and thus the 
distribution effectively corresponds to the Poisson dis- 
tribution. It is also important to note that log-normal 
distribution has an infinite maximum, unlike the log- 
binomial. To see this, recall that r(p) = —d p (p — 1), d p = 
D — D p , where D p are so-called generalized dimensions 
[22]. Therefore, 7 = d^ = lim^oo {-r(p)/(p - 1)}. For 
the log-normal process, d p = (p/2)p and 7 — ► 00. 

If we take the random [3 model "for real" , that is, con- 
sider the Poisson distribution as an approximation to the 
binomial, as in (11), then we are dealing with (6) with 
7/3 given from our measurements (so that /?2 is defined 
according to (7)). And the second parameter, the co- 
dimension Cp from (8), can be found from the best fit 
with the experimental data [7] , [8] - using computer rou- 
tines analogous to those we used for the Poisson distribu- 
tion (see above). The resulting Cp = 1.58, and we obtain 

for which 



„(bi) 



^ ESS) ) 2 ) = 0.0064, 



exponents Q 

quite satisfactory. Indeed, the corresponding r(p/3) de- 
picted in Fig. 3 is indistinguishable from other approxi- 
mations which collapse to the experimental data. 

In conclusion, one of the predictions of SL theory about 
the scaling of maxima ~ r -7 is experimentally confirmed. 
This makes it possible to make a better estimate of the 
intense structures geometry in fully developed turbu- 
lence. The PDF's of the exponents of the dissipation 



field are compared with the log-Poisson distribution to 
show a good agreement with the theory. The log-Poisson 
statistic can be considered as a limiting case for the log- 
binomial distribution appearing in random /3-model. We 
estimated the parameters of the log-binomial distribution 
with 7 found in our measurements, and to fit the expo- 
nents for the structure functions found elsewhere. We 
conclude that the estimated Hausdorff co-dimension of 
the most intense structures is less than 2. 

I thank K. R. Sreenivasan and B. Dhruva for sharing 
with me the data of atmospheric turbulence. I appreciate 
numerous comments made by S. Boldyrev, Z. Mikic, and 
R. Rosner. 



[1] A.N. Kolmogorov, C.R. Acad. Sci. U.S.S.R. 30, 301 
(1941). 

[2] A.N. Kolmogorov, J. Fluid Mech., 13, 82 (1962) 
[3] Z.S. She and E. Leveque, Phys. Rev. Lett. 72, 336 (1994). 
[4] B. Dubrulle, Phys. Rev. Lett. 73, 959 (1994). 
[5] Z.S. She and E.C. Waymire, Phys. Rev. Lett. 74, 262 
(1995). 

[6] B. Jouault, M. Greiner, and P. Lipa, Physica D, 136, 125 
(2000). 

[7] R. Bcnzi, S. Ciliberto, C. Baudet, R. Tripiccione, F. Ma- 
saioli, and S. Succi, Phys. Rev. E 48, R29 (1993). 

[8] R. Benzi, S. Ciliberto, C. Baudet, and Ruiz Chavarria, 
Physica D 80, 385 (1995). 

[9] R. Arratia, L. Goldstein, and L. Gordon, Stat. Sci. 5, No. 
4, 403 (1990). 

[10] K.R. Sreenivasan and R.M. Everson (unpublished). 
[11] S.I. Vainshtein, K.R. Sreenivasan, R.T. Pierrehumbert, 

V. Kashyap, A. Juneja, Phys. Rev. E, 50, 1823 (1994), 

see the end of Sec. VI. 
[12] R. Benzi, G. Paladin, G. Parisi, and A. Vulpiani, J. Phys. 

A, 17, 3521 (1984). 
[13] G.R. Ruetsch, and M.R. Maxey, Phys. fluids A 3 1587 

(1991); J. Jimenez, A. A Wray, P.G. Saffman, and R.S. 

Pogallo, J. Fluid Mech., 255, 65 (1993); H.K. Moffatt, S. 

Kida, and K. Ohkitani, J. Fluid Mech., 259, 241 (1994); 

T. Passot, H. Politano, P.L. Sulem, JR. Angilella, and 

M. Maneguzzi, J. Fluid Mech., 282, 313 (1955). 
[14] B.B. Mandelbrot, Lecture Notes in Mathematics, Vol 565 

page 121, Ed. R. Teman, Springer- Verlag (1976). 
[15] J. Clece, M. Greiner, and K.R. Sreenivasan, Europhysics 

Letters, 61, 756 (2003). 
[16] S. Boldyrev, A. Norlund, and P. Padoan, Astrophys. J. 

573, 678 (2002). 
[17] S. Boldyrev, Astrophys. J. 569, 841 (2002). 
[18] C. Meneveau and K.R. Sreenivasan, Phys. Rev. Lett. 59, 

1424 (1987). 

[19] U. Frisch, P-L. Sulem, and Nelkin, J. Fluid Mech. 87, 
719 (1978). 

[20] A.S. Monin and A.M. Yaglom, Statistical Fluid Mechan- 
ics, Vol. 2 (MIT Press, Cambridge, Mass, 1971) 



4 



[21] U. Frisch, Turbulence, (Cambridge University Press, 
!995). 

[22] H. G. E. Hentschel and I. Procaccia, Physica D 8, 435 
(1983). 



5 



